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We numerically investigate the transition of the static quark-antiquark string into a static-light meson-antimeson 
system. Improving noise reduction techniques, we are able to resolve the signature of string breaking dynamics 
for n/ =2 lattice QCD at zero temperature. We discuss the lattice techniques used and present results on energy 
levels and mixing angle of the static two-state system. We visualize the action density distribution in the region 
of string breaking as a function of the static colour source-antisource separation. The results can be related to 
properties of quarkonium systems. 


1. INTRODUCTION 

The breaking of the colour-electric string be¬ 
tween two static sources is a prime example of 
a strong decay in QGD p. Recently, we re¬ 
ported on an investigation of this two state sys¬ 
tem PEI!, with a wave function \Q) created by 
a QQ operator and a wave function \B) created 
by a four-quark BB operator, where B = Qq. Q 
denotes a static source and q is a light quark. 

We determined the energy levels Ei{r) — 2m b 
and E2{r) — 2m b of the two physical eigenstates 
|1) and 12) which we decomposed into the com¬ 
ponents, 

|1) = cos^ IQ) + sin6> |5) (1) 

|2) = — sin6> IQ) + cos |5). (2) 

We characterise string breaking by the distance 
scale Vc at which /S.E = E 2 —E 1 is minimized and 
by the energy gap AEc = AE{rc). While these 
energy levels and Vc are first principles QGD pre¬ 
dictions, the mixing angle 0 is (slightly) model de¬ 
pendent: within each (Fock) sector there are fur¬ 
ther radial and gluonic excitations and we trun¬ 
cated the basis after the four quark operator. 

In order to obtain dynamical information on 
the string breaking mechanism, we are study¬ 

* Based on talks presented by G.S. Bali and K. Schilling. 


ing the spatial energy and action density distri¬ 
butions within the two state system. In doing 
so one can address questions about the locali¬ 
sation of the light qq pair that is created when 
r is increased beyond Tc- The energy density 
will decrease fastest in those places where qq cre¬ 
ation is most likely. Perturbation theory suggests 
that light pair creation close to one of the static 
sources is favoured by the Goulomb energy gain 
while aesthetic arguments might suggest a sym¬ 
metric situation with qq dominantly being created 
near the centre. First results on this investigation 
are presented here. 

2. PARAMETERS AND METHODS 

We use Uf = 2 Wilson fermions at a quark 
mass slightly smaller than the physical strange 
quark and a lattice spacing a = 0.166(2) ro ~ 
0.083(1) fm and find jSj, Vc = 2.5(3) ro ~ 
1.25(1) fm and AEc ~ 51(3) MeV, where the er¬ 
rors do not include the phenomenological uncer¬ 
tainty of assigning a physical scale to ro ~ 0.5 fm. 
Using data on the QQ potential and the static- 
light meson mass m b , obtained at different quark 
masses, we determine the real world estimate, 
Vc = 1.13(10) (10) fm, where the errors reflect all 
systematics. An extrapolation of AEc however 


1 



2 


is impossible, without additional simulations at 
lighter quark masses. 

These results became possible by combining a 
variety of improvement techniques: the neces¬ 
sary all-to-all light quark propagators were calcu¬ 
lated from the lowest eigenmodes of the Wilson- 
Dirac operator, multiplied by 75 , after a vari¬ 
ance reduced stochastic estimator correction step. 
The signal was improved by employing a fat link 
static action. Many off-axis distances were imple¬ 
mented to allow for a fine spatial resolution of the 
string breaking region. For details see Ref. (Sj. 
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Figure 1 . Static-light effective masses, am^ = 
\r].[CB{t)/CB{t + a)], where CB{t) is the static- 
light mesonic correlation function. 
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Great care was also spent on optimizing the 
overlap with the ground states, within the \Q) 
and \B) sectors, using combinations of APE |5lbj 
and Wuppertal [3 smearing. Our APE smearing 
consists of the replacement, 

uif^^=Psuis){aUi^^+^^::l)- (3) 

Psu{ 3 ) denotes a projection operator, back onto 
the SU (3) group, and 


TjMjAn) rr(n)t 


( 4 ) 


i G {l,2,3},j G {±1,±2,±3}, is the spatial 

(n) 

staple-sum, surrounding /. We choose a = 2.0 

and define our APE smeared links, Ux,i = 

Eor the projection we use. 


PSU{3) (^) 
A' 


A' det(A')“^/^ 
^ e f/(3). 




e SU{3), 

( 5 ) 


The spatial transporters within our QQ states 
are products of APE smeared links, taken along 
the shortest lattice distance between the two 
sources. The APE smeared links are also em¬ 
ployed for the parallel transport within the Wup¬ 
pertal smearing of light quark sources, used to 
improve the static-light meson operators: 


4”+'^ = 


1 


1 + 6(5 


4 ”^+ <5 


An) 


±3 _ 

j=±l 


. ( 6 ) 


We set (5 = 4 and take the linear combination 
0 ( 20 ) _ 5.53230(^0) + 7 . 26040 ^^ 0 ) our smear¬ 
ing function. Best results are obtained by using 
smeared-local quark propagators. The quality of 
the overlap with the static-light mesonic ground 
state is visualized in the effective mass plot Eig- 
ure^ Note that we display the effective masses 
up to physical distances t > 1.2 fm. 

We were able to achieve values of jagp = 
K^-q|Q )|2 = 0.62(2) and \aB? = = 

0.96(1) at r Tc for the overlaps of our test wave 
functions |Tx) with the respective states on the 
right hand sides of Eqs. o and ©■ We display 
these results in Eigure [3 The almost optimal 
value ja^l ~ 1 was essential to allow Ei, E 2 and 
0 to be fitted from correlation matrix data. 


C{t) 


( C!QQ{t) Cgsit) \ 
V CBQ{t) CsBit) J 




obtained at moderate Euclidean times; t > 2a 
for Cbb and t > 4a for the remaining 2 matrix 
elements Cqq and Cqb = Cbq at r Tc- 
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Figure 2 . The overlaps of our test wave functions 
with the IQ) and \B) Fock states. String breaking 
takes place around Vc ~ 15 a. 


We follow Ref. (Hj and define action and energy 
density distributions, 

O'n(x) = ^ [^„(x) + H„(x)], (8) 

en(x) = ^ [Sn{x) - H„(x)] , (9) 

where 

An{x) = {n\A^(yi)\n) - {A^) (10) 

{$„(^)|A^(x,t/2)|$„(0)) , 

are our approximations to the creation opera¬ 
tors of the states |n), obtained from the diagonal- 
isation of C{t). We have suppressed the distance 
r from the above formulae and n = 1 denotes 
the ground state (dominantly QQ at r < rc) and 
n = 2 the excitation (dominantly BB at r < rc). 
Electric and magnetic fields are calculated from 
the plaquette, 

(^ + ^4) = ^ ^ , ( 11 ) 


Q _ _ 

B\x) = ^Y. [Ux,iAx)+U,-ai,ij ( 12 ) 

i=l 

B~Ux—aj,ij T Ux—a{i-\-j),ij'\ ’ 
where j = mod(i, 3) + 1 and, 

= ytr ' (13) 

We implement two different operators with the 
same continuum limits: in one case we identify 
Ux,^ with the link Ux,^ connecting x with x^afi. 
Additionally, we used smeared operators, 

Ux,^l = PsU{3) 

where 7 = 0.4 and is the sum over all six sta¬ 
ples enclosing in the three forward and in 

the three backward directions. The 7 -value was 
tuned to maximize the average plaquette, calcu¬ 
lated from the smeared links. For the un-smeared 
plaquette zq = 1 in Eq. while for smeared 
plaquettes zq = 1 0{as) is adjusted such that 

the vacuum expectation value of the average pla¬ 
quette remains unchanged. 



Eigure 3. The energy levels in physical units for 
Uf = 2. The bands represent the expected rif = 
2 + 1 scenario. 
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Figure 4. The mixing angle 0. 


The plaquette smearing enhances the sig¬ 
nal/noise ratio. Due to this smearing and the 
fat link static action used, the peaks of the dis¬ 
tributions around the source positions (that will 
diverge in the continuum limit) are less singular 
than in previous studies of SU{2) gauge theory 
at similar lattice spacings jHj. In the continuum 
limit the results from smeared and un-smeared 
plaquette probes will coincide, away from these 
self energy peaks. The draw back of plaquette 
smearing is that exact reflection positivity is vi¬ 
olated. However, our wave functions are suffi¬ 
ciently optimized to compensate for this. 

We insert the £’^(x, t) and H^(x, t) operators 
at position t/2 into the correlation matrix C(t), 
Eq. 0. For evenjodd t/a-values we average 
E^\B^ over the two adjacent time slices, respec¬ 
tively. Using the fitted ground state overlap ratio 
aqjaB and the mixing angle d as inputs, we cal¬ 
culate the action and energy density distributions 
Eqs. 0 and 0 in the limit of large t via Eq. m 
from the measured matrix elements. The distri¬ 
butions agree within errors within the time range 
3a < t < 6a. The results presented here are based 
on our t = Aa analysis. 

3. RESULTS 

To set the stage, we display the main results 
of Ref. ini in Figures IHI and 0 In the first figure 
we also speculate about the scenario in the real 


world with possible decays into as well as into 
BsBq. For our parameter settings and Uf = 2 
string breaking occurs at a distance Vc ~ 1.25 fm. 
In Figure 0 we show the mixing angle as a func¬ 
tion of the distance. The BB content of the 
ground state is given by sinO. Within our statis¬ 
tical errors 0 reaches 7r/4 at r = re- Remarkably, 
there is a significant four quark component in the 
ground state at r < Tc while for r > Vc the limit 
^ ^ 7r/2 is rapidly approached. 
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Figure 5. The transverse profile in the centre of 
the flux tube. 


Our basis vectors \Q) and \B) are no eigen¬ 
states of the system. The transition rate, g oc 
dCqs/dt, can be related to energy gap and mix¬ 
ing angle: g{r) ^ AE{r) sm{0) cos{0). This 
means that AEc ~ 2^(r). With a transition rate 
of only about 25 MeV in the string breaking re¬ 
gion and even smaller ^(r) -values at r > rc, a 
detection of the ground state contribution to the 
standard Wilson loop at r > Tc, where |1) con¬ 
tains little IQ) admixture, is hopeless: one would 
have to resolve the correlation function at times 
of order 1/g ^ S fm! For a complementary view 
on the problem, in the context of the breaking of 
higher representation strings in pure gauge theo¬ 
ries, see e.g. Ref. |S]. 
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Figure 6. Ground state action density distribution at r < re- 


Diagrammatically, g{r) can be represented as, 


lywvvi 



where Nc is the number of colours: the large 
Nc'> Uf expectation for the minimal energy gap 
is, AEc oc ^Jnf /Nc. Obviously, precision stud¬ 
ies at different Nc and n / are needed to establish 
the validity range of this prediction. The decay 
of the adjoint potential into two gluelumps, where 
AEc ^ 1 /^C 5 also fits nicely into this context. 

The quality of our density distribution data is 
depicted in Figure 0 for the ground state at a 
separation slightly smaller than Tc, as a function 


of the transverse distance x from the QQ axis. 
Due to cancellations between the magnetic and 
electric components the energy density is much 
smaller than the action density: for the compar¬ 
ison we have multiplied the energy density data 
by the arbitrary factor of 29. Note that the ra¬ 
tio cr/e will diverge like — In a A in the continuum 
limit. The differences between the shapes of the 
energy and action density distributions are not 
statistically significant. Below we only visualise 
the more precise action density results. 

We employ several off-axis separations. As¬ 
suming rotational symmetry about the interquark 
axis, each point is labelled by two coordinates. 
X denotes the distance from the QQ axis and y 
denotes the longitudinal distance from the cen- 
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Figure 7. Ground state action density distribution at r > Tc- 


tre point. We define an interpolating rectangular 
grid with perpendicular lattice spacing a and the 
longitudinal spacing slightly scaled, such that the 
static sources always he on integer grid coordi¬ 
nates. We then assign a quadratically interpo¬ 
lated value to each grid point z, obtained from 
points in the neighbourhood, |z — (x, ^)| < e = a. 
On the axis the data points are more sparse and 
we relax the condition to e = while for the 
singular peaks we maintain the un-interpolated 
values. 

In Figures Eland 13 we display ground state ac¬ 
tion density distributions for different distances 
around Vc- An mpeg animation has been pub¬ 
lished in Ref. 0 . The colour encodes the relative 
statistical errors and the lattice mesh represents 


our spatial resolution a. As already evident from 
the mixing angle 0{r) of Figure 0 above, string 
breaking takes place within a small region around 
Tc- All distributions are very similar to linear su¬ 
perpositions of string and broken-string states, 
with no non-trivial spatial dependence: string 
breaking appears to resemble an instantaneous 
process, without evidence of localisation of the 
qq pair creation. 

4. APPLICATION TO QUARKONIUM 
DECAY 

We wish to relate the static limit results to 
strong decay rates of quarkonia. In the non- 
relativistic limit of heavy quarks, potential “mod¬ 
els” provide us with the natural framework for 
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such studies. In fact at short distances, r <C 1/A, 
potential “models” can be derived as an effective 
field theory, potential NRQCD, from QCD pT7j . 
One can in principle add a BB sector, as well as 
transition terms between the two sectors, to the 
QQ pNRQCD Lagrangian. Strong decays would 
then be a straightforward non-perturbative gen¬ 
eralisation of the standard multipole treatment 
of radiative transitions in QED. Unfortunately, 
transitions such as T(45') ^ BB can hardly be 
classed as “short distance” physics. So, some 
modelling is required instead. 

The natural starting point again is a two chan¬ 
nel potential model which might still have some 
validity beyond the short distance regime: 


= Ei;{T) (16) 

with 

) p" + nr)- (17) 

tub / 

rriQ denotes the heavy quark mass and is the 
mass of a 5 meson. The wavefunction has two 
components. 



V.(r) = ( 

V V’bb(i’) 


and the potential is given by, 
E,{r) 


V{r) = 


0 

0 E2{r) 

Eqir) g{r) 
g\r) Eeir) 


O 


(18) 

(19) 

( 20 ) 


where we have normalized the zero point energy 
to 2mB and 


term g{r) and the BB continuum. In doing so, we 
assume that the interaction does not contain any 
spatial distribution but only depends on the dis¬ 
tance r. This instantaneous picture is supported 
by our action density measurements above. 



In Figure IHI we display our radial T(45') wave 
function u{r) = \/4^^'0qq 4 oo(^)* decay 

rate depends rather sensitively on the positions 
of the nodes. We obtain a preliminary result 
r[T(45') BB] ^ 5 MeV, which is about half 
the experimental value. This appears very rea¬ 
sonable, given the crudeness of the model and the 
fact that the gap AEc will increase with lighter, 
more realistic sea quark masses. We are studying 
the situation and systematics in more detail. 

5. CONCLUSIONS 


0=f “D (21) 

y —sm(J cos U J ^ ^ 

rotates our Fock basis {\Q),\B)} into the eigen- 
basis {|1), |2)}. We neglect B-B interactions and 
set Eb{t) = 0. We further adjust the difference 
^T( 4 S) “ 2mB to the experimental value which 
affects our phase space factor. We then follow 
Ref. cn and calculate the decay rate by multi¬ 
plying phase space with the overlap integral be¬ 
tween the T(45') wave function, the interaction 


We were able to resolve the string breaking 
problem in n/ = 2 QCD, at one value of the lat¬ 
tice spacing a~^ 2.37 GeV and of the sea quark 

mass, m < ms. It was also possible to study the 
dynamics of string breaking in detail and to re¬ 
solve spatial colour field distributions. The break¬ 
ing of the string appears to be an instantaneous 
process, with de-localized light quark pair cre¬ 
ation. While a direct lattice study of strong de¬ 
cay rates such as T(45') ^ BB or T(35') ^ DD 
is at present virtually impossible, investigations 




in the static limit can help constraining models. 
Studying the energy between pairs of static-light 
mesons can also be viewed as a milestone with re¬ 
spect to future calculations of AgAg forces, which 
are related to nucleon-nucleon interactions m 
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